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ORIGINAL ARTICLE 

The atypical ubiquitin ligase RNF3 1 stabilizes estrogen receptor a 
and modulates estrogen-stimulated breast cancer cell proliferation 

J Zhu\ C Zhao 1 , A Kharman-Biz 1 ' 2 , T Zhuang 1 , P Jonsson 3 , N Liang 1 , C Williams 3 , C-Y Lin 3 , Y Qiao 1 , K Zendehdel 2 , S Stromblad 1 , 
E Treuter 1 and K Dahlman-Wright 1 

Estrogen receptor a (ERa) is initially expressed in the majority of breast cancers and promotes estrogen-dependent cancer 
progression by regulating the transcription of genes linked to cell proliferation. ERa status is of clinical importance, as ERa-positive 
breast cancers can be successfully treated by adjuvant therapy with antiestrogens or aromatase inhibitors. Complications arise from 
the frequent development of drug resistance that might be caused by multiple alterations, including components of ERa signaling, 
during tumor progression and metastasis. Therefore, insights into the molecular mechanisms that control ERa expression and 
stability are of utmost importance to improve breast cancer diagnostics and therapeutics. Here we report that the atypical E3 
ubiquitin ligase RNF31 stabilizes ERa and facilitates ERa-stimulated proliferation in breast cancer cell lines. We show that depletion 
of RNF31 decreases the number of cells in the S phase and reduces the levels of ERa and its downstream target genes, including 
cyclin D1 and c-myc. Analysis of data from clinical samples confirms correlation between RNF31 expression and the expression of 
ERa target genes. Immunoprecipitation indicates that RNF31 associates with ERa and increases its stability and mono- 
ubiquitination, dependent on the ubiquitin ligase activity of RNF31. Our data suggest that association of RNF31 and ERa occurs 
mainly in the cytosol, consistent with the lack of RNF31 recruitment to ERa-occupied promoters. In conclusion, our study establishes 
a non-genomic mechanism by which RNF31 via stabilizing ERa levels controls the transcription of estrogen-dependent genes linked 
to breast cancer cell proliferation. 
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INTRODUCTION 

The relationship between estrogen signaling and breast cancer 
was revealed nearly 80 years ago. 1 The subsequent discovery and 
clinical application of antiestrogens and of aromatase inhibitors 
brought significant survival benefits to breast cancer patients. 2,3 
However, up to 40% of patients receiving adjuvant therapy 
eventually relapse, making resistance to endocrine therapy a 
significant clinical problem. 2,4 Thus, a detailed understanding of 
the underlying mechanisms and insight into new facets and 
components of estrogen signaling is critical in developing novel 
treatment strategies. 

Estrogen is known to exert its function by binding to estrogen 
receptor (ER) subtypes, ERa and ER(3, which belong to the nuclear 
receptor family of transcription factors. 5 ERa has a main role in 
breast cancer initiation and proliferation. Overexpression of ERa 
accelerates the G1-S phase transition, correlating with increased 
levels of oncogenic proteins, including cyclin D1 and c-myc. 6 These 
proteins promote cell cycle progression by decreasing the level of 
the cyclin-dependent kinase inhibitors p 2l Cip/WAF1 and p27 Kipl , the 
activity of which is essential for cell cycle progression. 7 Clinically, 
ERa levels in dysplastic breast epithelial cells correlate with the risk 
of breast cancer and response to endocrine therapy, and two- 
thirds of all breast cancers express high levels of ERa. 8,9 

There are a number of possible and confirmed mechanisms for 
the dysregulation of ERa that trigger inappropriate estrogen 
signaling and drug resistance in breast cancer. Some may relate to 



the transcriptional and epigenetic control of ERa expression, 10,11 
others to the control of ERa activity via ligands, growth factors and 
post-translational modifications. 12 " 14 However, our understanding 
of how ERa protein levels and stability are controlled in breast 
cancers remains largely unclear. As obvious candidates emerge 
components of the ubiquitin-proteasome system, including 
numerous E3 ubiquitin ligases, having been shown to have 
important roles in breast cancer by modulating estrogen signaling. 
For example, enzymes such as BRCA1, MDM2 and CHIP have been 
demonstrated to ubiquitinate ERa, suggesting that they directly 
trigger proteasomal degradation. 11,13,15,16 Despite these 
seemingly distinct mechanisms, a unifying model emerges 
wherein ERa turnover by the ubiquitin-proteasome system 
appears tightly linked to ER-mediated transcription. In fact, E3 
ligases appear to be recruited to ER-target genes, suggesting that 
ER ubiquitination and turnover occurs within short time frames at 
the chromatin and this tight coupling could be a prerequisite for 
subsequent rounds of transcription. 

That the regulation of ERa stability and turnover is complicated 
is evident from recent observations that other post-translational 
modifications such as phosphorylation, SUMOylation and acetyla- 
tion may crosstalk with ubiquitination. 11,17 Additionally, it is 
known that diverse ER-associated proteins can modulate ERa 
stability or even trigger proteasomal degradation and turnover, as 
demonstrated in the case of the oncogenic coactivator SRC3. 18 
Interestingly, ERa can be mono-ubiquitinated, which, unlike 
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poly-ubiquitination, does not trigger proteasomal degradation. 19 
However, the consequences of ERa mono-ubiquitination for its 
stability and function in breast cancer cells, as well as for the 
cellular factors that specifically trigger and recognize this type of 
modification, remain to be characterized. 

Here, we identify the atypical E3 ubiquitin ligase RNF31 (alias 
HOIP and ZIBRA) as one such candidate factor. RNF31 was initially 
cloned from breast cancer cells based on its elevated mRNA 
expression, 20 but its function in breast tumorigenesis and 
estrogen signaling has not yet been addressed. Our study 
characterizes a novel non-genomic mechanism that links RNF31 
to the control of ERa ubiquitination and stability and thereby to 
the transcriptional regulation of estrogen-dependent genes and 
breast cancer cell proliferation. 



RESULTS 

RNF31 modulates E2-stimulated proliferation of breast cancer cells 
To investigate the role of RNF31 in cell proliferation, we knocked 
down endogenous RNF31 in the ERa-positive breast cancer cell 
line MCF-7 using RNA interference (Figure 1a). We observed that 
RNF31 depletion reduced E2-stimulated cell growth (Figure 1b) 
and decreased E2-stimulated cell cycle progression (Figure 1c and 
Supplementary Figure S1 A). Importantly, RNF31 depletion reduced 
the number of cells in the S phase to the same extent as the 
depletion of ERa in E2-treated cells, which correlated with an 
increase in the number of cells in the G0-G1 phase. Over- 
expression of ERa partially reversed the reduced number of cells in 
the S phase following RNF31 depletion (Figure 1d and 
Supplementary Figure S1B). The actual level of reversion is 
difficult to interpret as the conditions for small interfering 
RNA (siRNA) transfection are different from those of plasmid 
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transfection and it is unlikely that all cells transfected with siRNA 
are also transfected with plasmid. Together these data suggest 
that ERa is a key mediator of the proliferative effects of RNF31 in 
this ER-responsive cellular system. However, also in the absence of 
E2 treatment (Figure 1c, vehicle), RNF31 depletion reduced the 
number of cells in the S phase, whereas ERa depletion had no 
effect under these conditions. Additionally, ERa depletion did not 
mimic RNF31 depletion with regard to the percentage of cells in 
the G0-G1 and G2-M phases of the cell cycle. 



RNF31 depletion decreases ERa levels and modulates expression 
of ERa target genes 

We then addressed the functional consequences of RNF31 
depletion on ERa signaling linked to proliferation control. Western 
blot analysis revealed that RNF31 depletion significantly down- 
regulated protein levels of ERa in MCF-7 cells (Figure 2a). 
Importantly, ERa downregulation was observed in a similar manner 
in MDA-MB175 and T47D cells (Supplementary Figure S2A) and 
could be achieved using two different individual siRNAs of the 
siRNA pool (Supplementary Figure S2B). To determine whether 
effects of RNF31 on ERa protein levels were correlated with effects 
on ERa transcriptional activity, we assayed ERa reporter gene 
activity following RNF31 depletion or overexpression. Figure 2b 
shows that RNF31 depletion leads to reduced ERa reporter gene 
expression, whereas RNF31 overexpression leads to increased 
reporter gene expression in the presence and absence of estrogen. 
RNF31 depletion also reduced the expression of endogenous ERa 
target genes such as ADORA 1, pS2 and Cydin Dl (Figure 2c). 
Consistent with this, chromatin immunoprecipitation analysis 
revealed decreased ERa binding to the promoter regions of target 
genes following RNF31 depletion (Figure 2d). Supplementary 
Figure 3A shows that inhibition of RNF31 does not affect the 



RNF31 
GAPDH 




siControl siRNF31 24h 48h 72h 

Time after transfection (h) 




Vehicle E2 Vehicle E2 Vehicle E2 G1 S G2-M 



Figure 1. RNF31 depletion inhibits cell proliferation and increases G1 arrest in MCF-7 cells, (a) MCF-7 cells were transfected with siRNF31 or 
siControl and knockdown efficacy was determined by western blot analysis for RNF31, using GAPDH as internal standard for control (Left 
panel), and qPCR (Right panel), (b) The WST-1 assay was used to determine the cellular metabolic activity at indicated time points after 
transfection. Cells are treated for indicated times with 10nivi E2 or vehicle. Experiments were done in triplicates. *P<0.05; **P<0.01 for 
siRNF31 E2 versus siControl E2. (c) RNF31 knockdown induces arrest in the G1 phase and inhibits the E2-mediated progression into the S 
phase. The effects of RNF31 knockdown were compared with the effects of ERa knockdown in MCF-7 cells. Cells are treated for 24 h with 1 0 nivi 
E2 or vehicle. The proportion of cells in each phase was measured by fluorescent-activated cell sorting (FACS). Experiments were done in 
triplicates. *P<0.05; **P<0.01; ***P< 0.001 for siRNF31 versus siControl. All values are mean ±s.d. (n = 3). (d) Overexpression of ERa partially 
reverses the reduced number of cells in the S phase following RNF31 depletion. MCF-7 cells were cultured in 10% fetal bovine serum (FBS). 
The proportion of cells in each phase was measured by FACS. Experiments were done in triplicates. *P<0.05; **P<0.01 for siRNF31 versus 
siControl; siRNF31 versus siRNF31 plus ERa overexpression. All values are mean±s.d. (n = 3). 
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Figure 2. RNF31 depletion decreases ERa protein levels and ERa signaling, (a) RNF31 depletion reduces ERa protein levels. MCF-7 cells were 
transfected with siRNF31 or siControl and treated with 1 0 n/vi E2 or vehicle for 72 h. ERa and RNF31 levels were determined by western blot 
analysis. GAPDH was used as internal control, (b) RNF31 depletion or overexpression affects ERa-dependent expression of an ERE-luciferase 
reporter gene. MCF-7 cells were transfected with siRNF31 or siControl or with plasmids expressing Myc-tagged RNF31 or Myc-tag vector alone 
or together with the ERE reporter plasmid. Subsequently, cells were treated with 1 0 nivi E2 or vehicle. Luciferase activity was measured 48 h 
after transfection. Shown are data from triplicate measurements. ***P< 0.001 for siRNF31 versus siControl and RNF31 overexpression versus 
control, (c) RNF31 depletion reduces the expression of endogenous ERa target genes. MCF-7 cells were transfected with siRNF31 or siControl. 
48 h after transfection, cells were treated with 1 0 nivi E2 or vehicle for 6 h. The expression levels of RNF3 7 and of endogenous ERa target genes 
{ADORA1, pS2, cyclinDI) were determined by qPCR from triplicate experiments. **P<0.01; ***P< 0.001 for siRNF31 versus siControl. (d) RNF31 
depletion decreases ERa recruitment to endogenous target gene promoters. MCF-7 cells were transfected with siRNF31 or siControl. Forty- 
eight hours post-transfection, cells were treated with 10 nivi E2 or vehicle for 30min and chromatin immunoprecipitation (ChIP) assays were 
performed with ERa antibody or rabbit immunoglobulin G (IgG) and quantified by qPCR. **P<0.01; ***P< 0.001 for siRNF31 versus siControl. 
(e) Heat map of ERa-regulated genes changed by RNF31 depletion in MCF-7 cells. P< 0.001 and fold change >2 was set as cutoff to derive 
regulated genes. All values are mean±s.d. (n = 3). 



endogenous expression of JUND and GAPDH, which were used as 
negative controls. Furthermore, Supplementary Figure S3B shows 
that inhibition of RNF31 effects ERa and nuclear factor-KB (NF-kB) 
signaling but not Liver X Receptor signaling in luciferase assays. 
Additionally, Supplementary Figure S3C shows that the lack of 
effect on Liver X Receptor signaling is independent of the presence 
of ligand. Thus, the effect of RNF31 on cell signaling shows 
pathway selectivity. The effect on the NF-kB pathway is not 
surprising considering the established role of RNF31 in modulating 
this pathway. 21-23 Consistent with the well-known regulation by 
ERa of its own expression, 24 RNF31 depletion downregulated the 
expression of ERa mRNA (Supplementary Figure S3D) and the 
binding of ERa to the known ERa-binding site in the ERa promoter 
(Figure 2d). Global gene expression analysis followed by sub- 
network enrichment analysis revealed significant regulation of ERa 
signaling pathways by RNF31 (Table 1). In line with this, RNF31 
affects a large number of ERa target genes, both those that 
have been shown to be upregulated and those that have been 
shown to be downregulated, in breast cancer cells (Figure 2e). 
Thus, RNF31 constitutes a regulator of general ERa signaling and 
its target genes. 



Table 1. Top 10 signaling pathways changed by RNF31 depletion in 
MCF-7 cells as determined by sub-network enrichment analysis 



Node 


Number of genes 


P-volue 


TNF 


119 


6.01 E - 08 


Sp1 


95 


3.49E - 07 


ERa 


40 


5.61 E- 07 


IL-ip 


79 


6.80E - 07 


TGF-P1 


100 


6.36E - 06 


TIMP3 


8 


9.35E-06 


EGF 


58 


1.43E-05 


MAPK3 


38 


1.69E-05 


GR 


35 


3.86E - 05 


SMAD7 


15 


5.33E-05 



Abbreviations: ERoc, estrogen receptor oc; EGF, epidermal growth factor; 
IL-1J3, interleukin-1 (3; MAPK3, mitogen-activated protein kinase 3; TGF-J31, 
transforming growth factor-pi; TNF, tumor necrosis factor; IMP3, tissue 
inhibitor of metalloproteinases-3. 
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Figure 3. RNF31 is highly expressed and is correlated to ERa target genes in tumour samples, (a) RNF31 is highly expressed in breast tumors. 
The expression of RNF31 was determined for 72 breast tumors and 37 adjacent breast tissues using qPCR. The Student's f-test was used for 
comparison. Whiskers, 10th and 90th percentile; box boundaries, 75th and 25th percentile; line within box, median. Dots above the boxes 
show sample maximum values, asterisk (*) indicates outliers, (b) Correlation between expression of RNF31 and expression of the classical 
positively regulated ERa target genes GREB1 and pS2 (TFF1). Expression levels are derived from the TCGA mRNA database, test statistics from 
Pearson product-moment correlation, (c) Heat map showing correlation between RNF31 mRNA levels and mRNA levels of known ERa target 
genes. Included ERa target genes correspond to genes displaying significant correlation with RNF31 in Supplementary Table S2. RNA 
expression data are from the TCGA RNA-sequencing database. 



RNF31 is highly expressed and is correlated to ERa target genes in 
tumor samples 

To begin to explore the clinical relevance of the effect of RNF31 on 
proliferation and estrogen signaling, we analyzed primary breast 
cancer samples and adjacent tissues for the expression of RNF31 
mRNA. We observed high levels of RNF31 expression in tumor 
tissue compared with adjacent tissues (Figure 3a). Next, we tested 
known ERa target genes that were also identified as being 
regulated by RNF31 in MCF-7 cells, for correlation with RNF31 in 
publically available gene expression profiling data from >2000 
breast cancer patients in the TCGA RNA-sequencing 25 and 
KMplot 26 databases. Importantly, RNF31 expression correlates 
with expression of about 70% of the known ERa target genes, in at 
least one of the clinical gene profiling data sets, identified as being 
regulated by RNF31 in MCF-7 cells (Table 2). This is exemplified in 
Figure 3b that shows that the expression of RNF31 positively 
correlates with the expression of classical ERa target genes, such 
as TFF1 (pS2) and GREB1 in clinical breast cancer samples. If 
tumors were separated according to high or low RNF31 
expression, a significant correlation to the differential expression 
of ERa target genes was observed, as further shown in Figure 3c. 

RNF31 associates with ERa and increases its stability 
Further support for the functional cooperation of RNF31 with 
ERa was obtained by co-immunoprecipitation (co-IP) of the 
overexpressed proteins from HEK-293 cells (Supplementary 
Figure S4A) and of the endogenous proteins from MCF-7 cells 
(Figure 4a and Supplementary Figures S4B and C). Additionally, 
overexpression of RNF31 upregulated ERa protein levels in MCF-7 
cells under steady-state conditions (Figure 4b). Upon inhibition of 
protein synthesis by cycloheximide, RNF31 depletion accelerated 
ERa degradation and this effect was observed both in the 



presence and absence of estrogen (Figures 4c and d). It is well 
established that ERa regulates its own expression in MCF-7 cells, 24 
making it difficult to distinguish direct effects of RNF31 on ERa 
protein and mRNA levels in this cell line. We therefore performed 
assays in HEK-293 cells, in which RNF31 but not ERa is expressed. 
In the presence of the proteasome inhibitor MG132, ERa was 
stabilized also in the absence of RNF31 and overexpression of 
RNF31 did not further increase the stability of ERa (Figure 4e). In 
summary, these data suggest that RNF31 increases ERa protein 
stability via inhibition of its proteasome-dependent degradation. 

The RBR domain of RNF31 is required for the association with, and 
stabilizing effect on, ERa 

To delineate the domains of RNF31 that are important for the 
association with ERa, wild-type RNF31 and deletion variants 
thereof (Figure 5a) were expressed together with ERa in HEK-293 
cells. Co-IP revealed a requirement of the C-terminal ring between 
ring fingers (RBR) domain for association with ERa (Figure 5b). We 
also overexpressed the same RNF31 variants in MCF-7 cells to 
determine the effect on endogenous ERa. In line with the notion 
that the RBR domain is required for RNF31 to associate with ERa, it 
is also required for the regulation of ERa levels (Figure 5c). In 
contrast, the ubiquitin-associated domain seems to be dispen- 
sable for associations with, and stabilization of, ERa. Finally, 
estrogen response element (ERE)-luciferase assays show that the 
RBR region of RNF31 is required to increase transcription- 
dependent ERa signaling in MCF-7 cells (Figure 5d). 

RNF31, via its RBR domain, triggers ERa mono-ubiquitination 

The requirement of the RNF31 RBR domain, implicated in 

catalyzing atypical ubiquitination of substrates, led us to 

investigate the consequences of RNF31 depletion or 
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Table 2. RNF31 expression correlates with expression of known ERa 
target genes 



Symbol 


TCGA 


Symbol 


KMplot.com 




All samples 




All samples 




r 2 


P-value 




r 2 


P-value 


ALb 


00.03 


00.37 


ALb 


00.25 


00.00 


ADA 


00.03 


00.31 


ADA 


00.06 


00.01 


AUUKA I 


00.04 


00.19 


AUUKA I 


00.11 


00.00 


A D/~\ A 1 

ArUA I 


00.17 


00.00 


A D/~\ A 1 

ArUA I 


00.08 


00.00 


AREG 


00.07 


00.05 


AREG 


NA 


NA 


A I IDI/ A 

AUKKA 


00.08 


00.01 


A I IDI/ A 

AUKKA 


-0.09 


01.00 


BTG2 


00.06 


00.08 


BTG2 


00.14 


00.00 


CCNG2 


00.02 


00.53 


CCNG2 


-0.10 


01.00 


/"CC1 D 

Lbrl K 


00.01 


0.05 


/"CC1 D 

Lbrl K 


00.21 


00.00 


rvni A 1 

LYrl Al 


00.02 


0.04 


f VD1 A 1 

LYrl Al 


00.22 


00.00 


CCD1 

tbK I 


00.18 


00.00 


CCD1 

tbK I 


00.13 


00.00 


LJdZ 


00.00 


0.06 


r~ id") 

LJbz 


NA 


NA 


GREB1 


00.21 


00.00 


GREB1 


00.11 


00.00 


i_i ^ /l/^v 1 
HIVIUA I 


00.01 


0.05 


HIVIUA I 


00.01 


00.37 


irrnn/i 


00.14 


00.00 


irrnn/i 

ILrbr4 


00.20 


00.00 


ixy- kiy a 1 

KLNIVlAl 


00.07 


00.03 


IX/— Ml/1 A 1 

KLNIVlAl 


00.07 


00.00 


rDZA I 


00.18 


00.00 


rUZK I 


NA 


NA 


bLKJ 


00.03 


00.34 


bbKJ 


NA 


NA 


yrri 

1 rr 1 


00.24 


00.00 


1 rr 1 


00.06 


00.00 


1 oroj 


00.10 


00.00 


1 Lrbo 


00.17 


00.00 


ACE2 


-0.17 


00.00 


ACE2 


-0.01 


0.04 


ATDD 

ALrr 


-0.01 


0.04 


Arnn 
ALrr 


00.18 


00.00 


LAVz 


-0.30 


00.00 


(~ A\/-) 

LAVz 


-0.08 


01.00 


CCNE1 


-0.10 


00.00 


CCNE1 


-0.09 


01.00 


LUoD 


-0.16 


00.00 


LDoO 


-0.06 


01.00 


LUKN I A 


-0.02 


00.57 


LUKN I A 


00.14 


00.00 


CXCL1 


-0.12 


00.00 


CXCL1 


-0.02 


0.05 


LXLL1 z 


-0.09 


00.01 


LXLL1 z 


00.00 


00.45 


LXLK4 


-0.02 


00.58 


LXLK4 


-0.10 


01.00 


rvrD7 
LXLK/ 


-0.17 


00.00 


LXLK/ 


-0.03 


0.06 


CYP19A1 


-0.07 


00.02 


CYP19A1 


00.16 


00.00 


CYP1B1 


-0.18 


00.00 


CYP1B1 


-0.07 


01.00 


EGFR 


-0.21 


00.00 


EGFR 


00.02 


00.17 


EGR3 


-0.01 


0.06 


EGR3 


00.10 


00.00 


NR5A2 


-0.14 


00.00 


NR5A2 


00.17 


00.00 


PPARA 


-0.20 


00.00 


PPARA 


00.13 


00.00 


REL 


-0.07 


00.02 


REL 


-0.01 


0.04 


S100A8 


-0.15 


00.00 


S100A8 


-0.10 


01.00 


TRPV6 


-0.09 


00.00 


TRPV6 


-0.03 


0.06 


WISP2 


-0.03 


00.39 


WISP2 


00.06 


00.01 



Abbreviations: ERa, estrogen receptor a; NA, not applicable. Yellow: 
P-value < 0.05; red: positive correlation; green: negative correlation. 



association with, and stabilizing effect on, ERa, the presence of 
mono-ubiquitinated ERa was dependent on the RBR domain of 
RNF31 (Figure 6e). To further support the functional link between 
mono-ubiquitination, stability and signaling by ERa, we used the 
ubiquitin ligase-deficient RNF31 mutant R1/2M, in which cysteine 
residues required for the transfer of ubiquitin to substrates have 
been mutated. 27 We overexpressed the RNF31 R1/2M in MCF-7 
cells to determine the effect on endogenous ERa and ERE- 
luciferase activity. Figures 7a and b show that R1/2M failed to 
increase endogenous ERa protein levels and ERE-luciferase 
activity. If co-transfected with ERa into HEK-293 cells, RNF31 
R1/2M failed to induce mono-ubiquitination of ERa (Figures 7c 
and d). Altogether, these data point at functional links between 
the RNF31 ubiquitin ligase domain, ERa mono-ubiquitination and 
ERa stabilization. 

RNF31 and ERa co-occur mainly in the cytosol 
RNF31 is known to be involved in diverse ubiquitination pathways 
both in the cytoplasm and in the nucleus, but its intracellular 
localization and action in relation to ER signaling in breast cancer 
cells remained unclear. Immunocytochemistry revealed that 
endogenous RNF31 is localized to the nuclear and cytosol 
compartments of MCF-7 cells (Figure 8a, for quantification see 
Supplementary Figure S5). Furthermore, endogenous ERa and 
RNF31 were found to be prevalent in the cytosol both in the 
presence and absence of E2 in MCF-7 cells (Figure 8a), and a similar 
protein distribution was seen with ectopically expressed proteins 
in HEK-293 cells (Supplementary Figure S6). In MCF-7 cells, 
treatment with E2 redistributed a significant fraction of ERa to 
the nucleus, whereas the RNF31 distribution remained unaffected. 
Consistent with the immunocytochemistry results, co-IP assays 
revealed that the association of RNF31 and ERa occurs mainly in 
the cytosolic fraction (Figure 8b). Furthermore, chromatin immu- 
noprecipitation assays failed to detect the recruitment of RNF31 to 
ER-responsive genes in MCF-7 cells (Supplementary Figure S7). 



DISCUSSION 

Here we report that the atyical ubiquitin ligase RNF31 associates 
with, and stabilizes, ERa in breast cancer cells, leading to increased 
estrogen signaling and cellular proliferation. Importantly, we also 
observe correlations between RNF31 and ERa target gene 
expression in breast cancer samples. Our study suggests a novel 
non-genomic mechanism of ER stability control (see model 
Figure 8c). On the basis of these data we propose that selective 
modulation of RNF31 expression and/or functional cooperation 
with ERa could be a strategy to inhibit proliferation in ERa-positive 
breast cancers. 



overexpression on the ubiquitination status of ERa in HEK-293 
cells and in MCF-7 cells. Western blot analysis of ERa in lysates 
from HEK-293 cells overexpressing RNF31 and ERa revealed the 
presence of potentially mono-ubiquitinated ERa species 
(Figure 6a). Importantly, depletion of endogenous RNF31 lead to 
reduced levels of mono-ubiquitinated ERa in MCF-7 cells 
(Figure 6b). The higher-molecular weight ERa species were 
reduced upon deletion of the RBR domain, whereas no change 
was seen upon deletion of the ubiquitin-associated domain 
(Figure 6c). More direct evidence that ERa was specifically subject 
to modification by ubiquitination was obtained by immunopre- 
cipitation-based ubiquitination assays. Ubiquitinated species were 
enriched using a ubiquitin antibody from HEK-293 cell extracts 
overexpressing ERa together with RNF31 variants and analyzed by 
western blots for the presence of ERa in the precipitates 
(Figure 6d). The results confirmed that RNF31 overexpression 
increased the amounts of mono-ubiquitinated ERa. Consistent 
with the important role of the RBR domain of RNF31 for its 



RNF31 — a mediator of crosstalk between proliferative signaling 
pathways in breast cancer cells 

Our analysis of MCF-7 cell proliferation, in conjunction with siRNA- 
mediated depletion of RNF31 or ERa, revealed a requirement of 
RNF31 for breast cancer cell growth and specifically for E2- 
stimulated cell cycle progression. This suggests that ERa is one of 
the key mediators of the proliferative effects of RNF31 in ERa- 
positive breast cancer cells. Indeed, unbiased pathway analysis of 
the RNF31 -dependent transcriptome, as determined by microarray 
expression analysis in MCF-7 cells, confirmed the functional link of 
RNF31 to the estrogen-signaling pathway. Importantly, a func- 
tional link between RNF31 and the estrogen-signaling pathway 
was also observed in breast cancer samples. In addition, the 
analysis revealed links to the NFkB pathway, and to the p53 
pathway, known for its key role in cell cycle control. In addition, a 
previous study reported that RNF31 participates in transforming 
growth factor (3 and Wnt/(3-catenin signaling pathways, 28 both of 
which play fundamental roles in breast cancer. These pathways 
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Figure 4. RNF31 associates with ERa and increases its stability, (a) Co-IP assays reveal associations between endogenous RNF31 and ERa in 
MCF-7 cells, (b) Over-expression of RNF31 increases endogenous ERa protein levels in MCF-7 cells. MCF -7 cells were transfected with 
plasmids expressing Myc-tagged RNF31 or the Myc-tag alone. Forty-eight hours after transfection, whole-protein extracts were prepared and 
subjected to western blot analysis. The levels of RNF31, ERa and the internal control GAPDH were determined by western blot analysis. The 
predicted molecular weights are indicated, (c) Depletion of RNF31 decreases endogenous ERa protein stability in MCF-7 cells. Cells were 
transfected with siRNF31 or siControl. After 72 h, cells were treated with 100jamol/l cycloheximide for different times as indicated before 
whole protein extraction. The levels of ERa and the internal control GAPDH were determined by western blot analysis. The experiments were 
done in triplicates. Image J was applied to quantify ERa and GAPDH density, (d) Depletion of RNF31 decreases endogenous ERa protein 
stability in MCF-7 cells in the presence of E2. Cells were transfected with siRNF31 or sicontrol. After 72 h, cells were treated with 1 0 nrvi E2 for 
24 h. Then cells were treated with 100jimol/l cycloheximide for different times as indicated before whole protein extraction. The levels of 
RNF31, ERa and the internal control GAPDH were determined by western blot analysis. The experiments were done with triplicates. Image J 
was applied to quantify the ERa and GAPDH density, (e) Overexpression of RNF31 does not further increase the stability of ERa in the presence 
of the proteasome inhibitor MG132. HEK-293 cell were transfected with plasmids expressing Myc-tagged RNF31 or the Myc-tag alone. Forty- 
eight hours post transfection, cells were treated with 100nmol/l MG132 or vehicle. Cells were harvested 6h after MGM132 treatment and 
whole protein extracts were prepared. The levels of RNF31, ERa and the internal control GAPDH were determined by western blot analysis. 



may contribute to the ERa/E2-independent effects of RNF31 
depletion, consistent with the observation that RNF31 is also 
expressed in ERa-negative cell lines (data not shown). 

Irrespective of the precise contribution of these additional 
RNF31 -regulated pathways, it is likely that they crosstalk with the 
estrogen signaling pathway in ERa-positive cancers, instead of 
acting independently. One intriguing candidate for such cross talk 
is the NFkB signaling pathway, implicated in both inflammation 
and oncogenesis in cancer cells. 29,30 A number of recent studies 
have demonstrated that RNF31, together with the linear ubiquitin 
chain assembly complex (LUBAC) complex, is a critical component 
of cytoplasmic NFkB activation cascades, in other cell types. 21-23 
However, a first direct evidence for both positive and negative 
crosstalk has been provided by studying tumor necrosis factor re- 
activation of the NFkB pathway, and its interference with ER 
signaling, in breast cancer cell lines. 29 Future studies of the 
individual facets of this crosstalk with particular emphasis on the 
RNF31 network should help to re-evaluate newly developed anti- 
inflammatory ER ligands, 29 along with attempts to target 
components of the RNF31/NF-kB pathway. 



The place of RNF31 within the atypical RBR family of E3 ubiquitin 
ligases 

Among the more than 700 putative E3 ubiquitin ligases in 
humans, members of the RING-in between rings (IBR)-RING (RBR) 
family have attracted recent attention due to their uncommon 



ubiquitination mechanism and their regulatory signaling 
functions. 31 Unlike classic E3s, including RING and HECT ligases, 
RBR ligases seem to preferentially catalyze mono-ubiquitination 
or linear poly-ubiquitination of substrates, which does not lead 
to substrate degradation and turnover. 32,33 The very recent 
clarification of the unique RBR ubiquitination mechanism in part 
resolves the initial paradox that RBR ligases often stabilize 
their substrates. Additionally, mechanisms of RBR-mediated 
ubiquitination are often substrate-dependent. 31,32 RBR ligases 
appear therefore well suited to fulfill diverse regulatory roles in 
intracellular signaling and transcription. This has perhaps best 
been demonstrated in case of the Parkin protein (PARK2), 
implicated in Parkinson disease, and in the case of PARC (CUL9) 
which is implicated in p53 tumor suppression. 34-37 

There is a growing interest in the linear ubiquitin chain 
assembly complex referred to as LUBAC, implicated in NFkB 
signaling, and thus potentially in multiple disorders, such as 
cancer and chronic inflammatory and metabolic diseases 21,22 
Most notably, RNF31 functions as the main E3 ligase of LUBAC, 
which contains additional subunits including RBCK1 and Sharpin. 
Both RNF31 and RBCK1 function independently from each other, 
and thus without the entire LUBAC complex, to modulate various 
nuclear receptor pathways. For example, we have previously 
demonstrated that RBCK1 inhibits cellular proliferation of ERa- 
positive breast cancer cells via several genomic/transcriptional 
mechanisms: RBCK1 regulates ERa gene transcription, and it acts 
as a transcriptional cofactor that co-occupies E2-responsive 
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Figure 5. The RNF31 RBR domain is required for association with and stabilization of ERa. (a) RNF31 domain structure and deletion mutants 
used in this study (full-length, AUBA and ARBR). ZnF_RBZ, putative zinc-finger ubiquitin-binding domain; UBA, ubiquitin-associated domain, 
mediates interaction with RBCK1/LUBAC; RING-IBR-RING (RBR), atypical E3 ubiquitin ligase domain, (b) RNF31 association with ERa requires the 
RBR domain. HEK-293 cells were transfected with ERa together with plasmids expressing Myc-tagged full-length RNF31, variants deleting the 
UBA or RBR domains, respectively, or the Myc-tag alone. Cells were harvested 24 h after transfection and whole-cell extracts were prepared for 
co-IP. The predicted molecular weights of the RNF31 derivatives and ERa are indicated, (c) The RBR domain is necessary for the RNF31- 
mediated increase of endogenous ERa protein levels. MCF-7 cells were transfected with plasmids expressing Myc-tagged RNF31 derivatives or 
the Myc-tag alone, as indicated. 48 h after transfection, whole-cell extracts were prepared and levels of ERa protein assayed by western blot 
analysis. The predicted molecular weights of RNF31 variants, ERa and the loading control GAPDH are indicated, (d) The RBR domain of RNF31 
is required to increase ERa signaling. MCF-7 cells were transfected with plasmids expressing Myc-tagged RNF31 derivatives, or the Myc-tag 
alone, as indicated, along with an ERE-luciferase reporter plasmid. Luciferase activity was measured 48 h after transfection and calculated from 
experiments performed in triplicates. Data are shown as mean ±s.d. (n = 3). ***p< 0.001 for wt-RNF31 versus RNF31 deletion domains/empty 
vector; NA, P>0.05 for wt-RNF31 versus RNF31 deletion domains. 



promoters. 38 The combined data of our previous and current 
studies suggest that RNF31 and RBCK1 affect ERa signaling by 
very distinct mechanisms. Unlike RBCK1, RNF31 cannot be directly 
recruited to ERa target promoters in MCF-7 cells, consistent with 
the dominant cytoplasmic localization of RNF31 that we found. 

RNF31 stabilizes ERa via a unique ubiquitination mechanism 
In contrast to RBCK1, RNF31 appears specifically linked to the 
control of ERa protein stability and mono-ubiquitination. This 
mechanism of action is consistent with a previously reported and 
similar role of RNF31 in stabilizing the steroidogenic orphan 
receptor DAX-1. 27 The obvious link between RNF31 and mono- 
ubiquitination explains why RNF31 action leads to substrate 
stabilization instead of degradation, although some mechanistic 
details remain unrevealed. On the basis of our experimental data, 
the following two mechanisms are likely to contribute: (1) RNF31 
mono-ubiquitinates ERa, either alone or in complex with yet 
another E3 ligase. Importantly, RNF31 does not act in complex 
with RBCK1/LUBAC, as deletion of the RNF31 ubiquitin-associated 
domain, necessary for RBCK1 interaction, does not abolish mono- 
ubiquitination and stabilization of ERa (this study), or of DAX-1. 27 
However, the enzymatically intact RBR-type ubiquitin ligase 



domain is necessary for these events, and its activity might be 
further modified in a substrate (that is, ERa)-dependent manner. 
(2) RNF31 stabilizes ERa by preventing poly-ubiquitination, which 
may involve direct binding to mono-ubiquitinated ERa. 
Intriguingly, RNF31 has several putative Zn-finger type ubiquitin- 
binding domains, whose functionality and specificity, due to the 
lack of consensus signatures, 27 remain to be determined. 



RNF31 — a candidate target for breast cancer diagnosis and 
therapy? 

What evidence supports the clinical relevance of our findings 
about the role of RNF31 in estrogen signaling? Indeed, RNF31 was 
originally discovered in screens for genes that were differentially 
expressed in breast cancer cells compared with normal breast 
epithelial cells 20 Consistent with this, using a limited number of 
samples, we found that RNF31 is more highly expressed in breast 
cancer tissues as compared with the surrounding normal tissues. 
In addition, consistent with the expression of RNF31 mRNA in ERa- 
positive breast cancer cell lines, RNF31 mRNA appears highly 
expressed in ERa-positive breast tumor tissue. Importantly, the 
notion that RNF31 is expressed in both ERa-positive and negative 
cancers is not contradictory to our model but rather points at the 
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Figure 6. RNF31 triggers ERa mono-ubiquitination. (a) Detection of a potentially mono-ubiquitinated form of ERa upon RNF31 overexpression. 
HEK-293 cells were transfected with ERa together with plasmids expressing Myc-tagged RNF31 or the Myc-tag alone. Forty-eight hours after 
transfection, whole-cell extracts were prepared and levels of ERa protein assayed by western blot analysis. The predicted molecular weights of 
RNF31, ERa, mono-ubiquitinated ERa and of the internal control GAPDH are indicated, (b) Detection of endogenous mono-ubiquitinated ERa 
upon RNF31 depletion. MCF-7 cells were transfected with siRNF31 or siControl. Forty-eight hours after transfection, whole-cell extracts were 
prepared and levels of ERa protein assayed by western blot analysis. The predicted molecular weights of RNF31, ERa, mono-ubiquitinated ERa 
and the internal control GAPDH are indicated, (c) Deletion of the RNF31 RBR domain abolishes the potentially mono-ubiquitinated form of 
ERa. HEK-293 cells were transfected with ERa together with plasmids expressing Myc-tagged full-length RNF31 derivatives or the Myc-tag 
alone. Forty-eight hours post transfection, cell extracts were prepared and ERa forms were detected by western blot analysis. The predicted 
molecular weights of RNF31, ERa, mono-ubiquitinated ERa and of the internal control GAPDH are indicated, (d) Direct evidence for ERa mono- 
ubiquitination. Immunoprecipitation of ubiquitinated proteins from MCF-7 cell extracts upon overexpression of RNF31. Ubiquitinated ERa 
species were detected by western blots using anti-ERa, identifying a prominent 75 kDa mono-ubiquitinated ER form, (e) ERa mono- 
ubiquitination requires the RNF31 RBR domain. Plasmids expressing Myc-tagged RNF31 derivatives were transfected into HEK-293 cells 
together with the ERa expression plasmid. Whole-cell extracts were subjected to immunoprecipitation of ERa and subsequently analyzed for 
ubiquitinated ERa forms by western blot analysis using anti-ubiquitin. The predicted molecular weight of mono-ubiquitinated ERa is 
indicated. 



convergence of distinct signaling pathways toward proliferation 
and tumorigenesis, as discussed above. 

Given our data, it will be interesting but technically challenging 
to determine whether high RNF31 protein levels in tumor samples 
also cause a stabilization of ERa protein levels. Optimally, cancer 
samples containing equal levels of ERa mRNA but different levels 
of RNF31 (mRNA and protein) should be quantitated for ERa 
protein. An additional complication arises from the fact that 
functional associations between RNF31 and ERa occur mainly in 
the cytoplasm, which will be more complicated to analyze by 
immunohistochemistry than the nuclear fraction, traditionally 
used to score for ERa positivity in clinical samples. These concerns 
apply to monoclonal antibody-based biochemical methods as 
well, as already noted in 1986. 39 An interesting question is 
therefore whether some tumors that were initially scored as ER- 
negative could be enriched in cytoplasmic receptors, and whether 
this enrichment is dependent on cytoplasmic RNF31. Also relevant 



is the recent observation that clinically used markers including 
ERa, PR and HER2 protein levels are unstable throughout tumor 
progression and possibly influenced by adjuvant therapy. 9 As 
alterations of both epigenetic/transcriptional (mRNA) and post- 
translational (protein) mechanisms are likely to contribute to the 
loss of ERa, post-translational ER regulators such as RNF31 have to 
be taken into consideration as well. In light of these concerns, the 
possible convergence of RNF31 and ER signaling in RNF31/ERa- 
positive breast cancers should be of outmost interest for the 
clinical diagnosis and treatment decisions and requires further 
clarification. 

In conclusion, although ERa has already been well documented 
to have a critical role in the etiology and progression of breast 
cancer, RNF31 emerges as a new component of ER signaling 
in ERa-positive breast cancers. To identify novel therapeutic 
strategies for treating ERa-positive breast cancers, modulation of 
ERa levels is one feasible approach to inhibit estrogen signaling 
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Figure 7. The ubiquitin ligase activity of RNF31 is required for ERa stabilization, signaling and mono-ubiquitination. (a) The mono- 
ubiquitination function of RNF31 is required for the RNF31 -mediated increase of endogenous ERa protein levels. MCF-7 cells were transfected 
with plasmids expressing Myc-tagged RNF31, Myc-tagged RNF31 R1/2M, in which cysteine residues responsible for the transfer of ubiquitin to 
substrates have been mutated or the Myc-tag alone, as indicated. Forty-eight hours after transfection, whole-cell extracts were prepared and 
levels of ERa protein assayed by western blot analysis. The predicted molecular weights of RNF31 variants, ERa and the loading control 
GAPDH are indicated, (b) The mono-ubiquitination function of RNF31 is required to increase ERa signaling. MCF-7 cells were transfected with 
plasmids expressing Myc-tagged RNF31, Myc-tagged RNF31 R1/2M or the Myc-tag alone, as indicated, along with an ERE-luciferase reporter 
plasmid. Luciferase activity was measured 48 h after transfection and calculated from experiments performed in triplicates. Data are shown as 
mean±s.d. (n = 3). ***P< 0.001 for wt-RNF31 versus RNF31 R1/2M or Myc-tag. RNF31 deletion domains/empty vector; NA, P>0.05 for 
wt-RNF31 versus RNF31 deletion domains, (c, d) Mutation of the RNF31 ubiquitin ligase domain abolishes mono-ubiquitination of ERa. 
(c) HEK-293 cells were transfected with ERa together with plasmids expressing Myc-tagged RNF31, RNF31 R1/2M or the Myc-tag alone. Forty- 
eight hours post transfection, cell extracts were prepared and ERa was detected by western blot analysis. The predicted molecular weights of 
RNF31, ERa, mono-ubiquitinated ERa and the internal control GAPDH are indicated, (d) HEK-293 cells were transfected with ERa together with 
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The predicted molecular weights of RNF31, IgG, ERa and mono-ubiquitinated ERa are indicated. 



and subsequently cellular proliferation. As our study identifies 
RNF31 as a modulator of ERa protein levels in breast cancer cells, 
RNF31 could constitute a valuable diagnostic tool and/or a drug 
target for ERa-positive breast cancers. 

MATERIALS AND METHODS 

Cell culture 

MCF-7 and HEK-293 cells were cultured in Dulbecco's Modified Eagle 
Medium (DMEM) (Invitrogen, Carlsbad, CA, USA) supplemented with 10% 
fetal bovine serum and 1% penicillin/streptomycin (Invitrogen) at 37 °C in a 
humidified atmosphere of 5% C0 2 in air. T47D and MDA-MB-1 75 cells were 
cultured in RPMI 1640 (Invitrogen) supplemented with 10% fetal bovine 
serum and 1% penicillin-streptomycin. 1 7p-estradiol (E2; Sigma-Aldrich, 
St Louis, MO, USA) was dissolved in ethanol. For proliferation assays, cells 
were cultured in fetal calf serum treated with dextran-coated charcoal and 
phenol red-free DMEM. 

Plasmids 

RNF31 (pcDNA-Myc-RNF31) constructs were described previously. 27 The 
pcDNA3-ERa plasmid, the ERE-TK-Luc reporter and the pRL-TK control have 
been described previously. 40 



siRNA transfection 

Cells were transfected with 50-100 nivi siRNA. RNF31 siRNAs correspond 
to the siGENOME SMARTpool RNF31 (Dharmacon D-0021419, Thermo 
Fisher Scientific, Waltham, MA, USA), consisting of four single 
siRNAs, which were used in previous studies. 27 ERa siRNA and control 
siRNA were Stealth Select siRNA (Invitrogen). INTERFERin transfection 
reagent (Polyplus Transfection, New York, NY, USA) was used according to 
the manufacturer's protocol. 



Quantitative PCR 

Quantitative PCR was performed as previously described. 41 Primer 
sequences for qPCR are provided in Supplementary Table SI. 

Western blotting 

Cells were lysed with NP-40 lysis buffer (Roche Diagnostics GmbH, 
Mannheim, Germany). Antibodies used were anti-GAPDH (ab8425), anti- 
RNF31 (ab46322) from Abeam, and anti-Myc (9E10), anti-ubiquitin (FL76), 
anti-HA antibody (F-7), anti-ERa (1D5 and HC-20) from Santa Cruz 
Biotechnology, Inc., Dallas, TX, USA. 
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Figure 8. Intracellular localization analysis and model of crosstalk between RNF31 and ERa signaling in breast cancer cells, (a) MCF-7 cells were 
treated with 10nM E2 or vehicle for 30min before fixation. Intracellular localization of RNF31 (red) and ERa (green) was determined by 
immunofluorescence staining. Nuclei (blue) were stained with 4',6-diamidino-2-phenylindole (DAPI). Shown are representative images. For 
quantitative analysis see Supplementary Figure S5. (b) Co-IP assay reveals the interaction between RNF31 and ERa in the cytoplasm. The 
Subcellular Protein Fractionation Kit (Thermo Scientific, 78840) was used for extraction of cytoplasmic and nuclear proteins from MCF-7 cells. 
Vinculin and Histone 3 were used to identify the quality of cytoplasmic and nuclear fractions, respectively (left panel). Co-IP assays were 
performed with RNF31 antibody for precipitation and ERa antibody for detection (right panel), (c) Hypothetical model for the functional 
interplay of RNF31 with ERa signaling in breast cancer cells. RNF31 associates with ERa predominantly in the cytoplasm and promotes mono- 
ubiquitination of ERa, thereby counteracting proteasome-mediated degradation. The so elevated ERa protein levels cause increased genomic 
ERa signaling, for example, transcription of E2-target genes linked to the proliferation of breast cancer cells. RNF31, as part of the LUBAC 
ubiquitin ligase complex, has an established independent function in NFkB signaling. To which extent the two signaling pathways also 
crosstalk remains to be investigated. 



Luciferase reporter assays 

Luciferase activity was measured using the Dual-Luciferase Reporter Assay 
(Promega, Mannheim, Germany). 

Quantification of cell viability 

The cell number was determined using the WST-1 cell proliferation reagent 
(Roche) as previously described. 38 

Fluorescent-activated cell sorting 

MCF-7 cells were transfected with siRNF31, siERa or siControl. After 24 h, 
the medium was changed into phenol red-free DMEM with 2% dextran- 
coated charcoal-treated fetal calf serum. After an additional 24 h, the cells 
were fixed in 70% ethanol for 30 min and stained by propidium iodide. The 
BD LSR II flow cytometer (BD Bioscience, Franklin Lakes, NJ, USA) was used 
to measure the flow fluorescence intensity. For the rescue experiment, 
MCF-7 cells were transfected with siRNF31 or siControl. After 24 h, the cells 
were transfected with 2|ig Flag-ERa vector or empty Flag vector and 
cultured in DMEM medium with 10% fetal calf serum. 

Chromatin immunoprecipitation assays 

Chromatin immunoprecipitation assays were performed essentially as 
previously described. 41 For detail see Supplementary Data. 

Co-immunoprecipitation 

Co-IP was performed essentially as previously described. 41 The Subcellular 
Protein Fractionation Kit (Thermo Scientific, Waltham, MA, USA, 78840) was 
used for extraction of cytoplasmic and nuclear proteins. 



Protein stability assay 

MCF-7 cells were transfected with lOOnivi siRNF31 or siControl. 24 h post- 
transfection, cells were treated with E2 or vehicle for 24 h. The cells were 
then treated with cycloheximide (100(im) or MG132 (100um) for 48 h. 
Samples were harvested and analyzed by western blotting. 

Ubiquitination assays 

(A) To detect modified ERa from cell extracts, HEK-293 cells were 
transfected with 4[ig pcDNA3-ERa together with pCDNA3-Myc-RNF31 
including deletion variants. Forty-eight hours post-transfection, cells were 
treated with IOuM MG132/vehicle for 8h, directly lysed and modified/ 
unmodified ERa was detected by western blotting. (B) To directly enrich 
ubiquitinated ERa from cell extracts, HEK-293 cells were transfected with 
4|ig pcDNA3-ERa together with 4ug pCDNA3-Myc-RNF31 or deletion 
variants. Forty-eight hours post-transfection, total protein was extracted 
and pre-cleared by 20|il protein A slurry and 1 |ig rabbit immunoglobulin 
G for 2h at 4°C. The supernatant was collected and immunoprecipitated 
with 20jil ubiquitin beads (CalBiochem, no. 662200). Western blot with 
rabbit anti-ERa antibody was performed to detect ubiquitinated ERa. 

Immunofluorescence assay 

A rabbit anti-RNF31 polyclonal antibody and mouse anti-ERa monoclonal 
antibodies were used. For detail see Supplementary Data. 

Microarray analysis 

The Agilent SurePrint 8 x 60K arrays were used for global gene expression 
analysis at the core facility for bioinformatics and expression analysis (BEA, 
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www.bea.ki.se), Karolinska Instituted Sweden, according to standard 
protocols. Data are deposited in the Gene Expression Omnibus (GEO) 
database (accession number GSE46010). We applied a filter of P< 0.001 for 
significantly modulated gene expression and at least a 2.0-fold change in 
mean differential expression. 

Analysis of human breast tumor samples 

Human samples were from the Iranian Tumor Bank, Cancer Institute of Iran. 
For detail see Supplementary Data. 

Analysis of gene expression in publicly available data sets 
Pathway analysis and identification of regulated ERa target genes from 
microarray data were carried out using Pathway Studio (Elsevier). Analyses 
of gene expression in breast cancer samples from TCGA 25 and 
KMplot.com 26 were carried out in the statistical environment R. 
Presented correlations are Pearson's product-moment correlations. 
Samples were separated into high/low expression of RNF31 by setting a 
cutoff at the median value. 

Statistics 

Student's f-test, Pearson correlation coefficient and Cox regression analysis 
were used for comparisons. A P-value of <0.05 was considered to be 
significant. 
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